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This work confirms that not only surface plasmons but many other kinds of electromagnetic 
eigenmodes should be considered in explaining the values of the transmittivity through a slab bearing 
a two-dimensional periodic corrugation. Specifically, the role of Brewster-Zennek modes appearing 
in metallic films exhibiting regions of weak positive dielectric constant. It is proposed that these 
modes play a significant role in the light transmission in a thin chromium film perforated with 
normal cylindrical holes, for appropriate lattice parameters. 

PACS numbers: 78.20.-e, 42.79.Dj, 42.25.Bs, 42.25.Gy 



I. INTRODUCTION 

The optical properties of thin metallic films contain- 
ing periodic arrays of subwavelength holes have been ac- 
tively studied in relation with the possibility of trans- 
mission engineering. Originally, the remarkable optical 
transmission found in these structures was pointed out 
by T.W. Ebbesen et al back in 1998i. After these early 
experimental observations, the role of the thin metal- 
lic film surface plasmons (SPs) was put forward in or- 
der to explain the peculiar wavelength dependence of the 
transmissionii2i2i^iS*S*2iSiS. As a consequence of this sug- 
gestion, it was generally believed that the typical trans- 
mission features could only be obtained with metallic 
filmsii^. However it is important to note that other ex- 
planations had been developed which could compete with 
the SPs niode&iSiiiiiSiiiii. For instance, it was sug- 
gested that these phenomena could also be described in 
terms of the short range diffraction of evanescent waves 
or in terms of dynamical diffraction cffccts^^. Another 
explanation suggests to emphasize the role of cavity res- 
onances taking place in the holes to explain the trans- 
mission enhancement^. 

Many questions remain to completely clarify the scat- 
tering processes involved in these experiments. Nev- 
ertheless, even if the exact role of SPs is not clearly 
assessecli*2iMiSiSiL2iS, most interpretations tend to ad- 
mit that SPs should play a key role in the observed 
transmission spectrum features. The observed transmis- 
sion exhibits a set of peaks and dips. Many authors 
suggested that transmission peaks are created by SPs 
resonance3iiSiiiiSiSiL& but others rather saw the SPs res- 
onances at the spectral dips^. In a recent worki^ we sug- 
gested to redefine the role of SPs in Ebbesen experiments 
in the required context of resonant Wood anomalies^^. In 
doing so, we have shown that the transmission spectrum 
could be better depicted as a series of Fano profilesi^. 



FIG. 1: A view of the system under study 



These recognizable lineshapes result from the interfer- 
ence of non resonant transfers with resonant transfers 
which involve the film SPs and evanescent diffraction or- 
ders. We could then point out that each transmission 
peak-dip pair is nothing else than a Fano profilei^. The 
appearance of an assymetric Fano lineshape does not nec- 
essarily locates the peak or the dip at the SP resonance!^. 
However, we have shown that the existence of the SPs is 
a condition for the presence of the Fano lineshapesi^. It 
can be noted that recent results by C. Genet et al confirm 
this description^**. 

As a further outcome of this work-'^^ it became clear this 
kind of transmission spectrum, with its Fano lineshapes, 
could also be obtained in a more general context and that 
the generic concept of eigenmodes could substitute that 
of SPs. These results implied two important new ideas. 
First, according to the contributions of resonant and non- 
resonant processes in the Fano profile, eigenmodes can be 
associated with wavelengths closer to the peaks or to the 
dips. Second, it is possible to obtain transmission curves 
similar to those of metal films, by substituting SPs with 
guided modes or other types of polaritons. Many exam- 
ples can be found, including highly refractive materials 
defining guided modes or ionic crystals in the restrahlen 
band defining phonon-polaritons. 

In a previous paper*'?, it was made clear that a strong 
hypothesis on the origin of the eigenmodes was not re- 
quired to account for the general features of the mem- 
brane transmission. As metal films were considered 
there, it seemed natural to involve SPs in the transmis- 
sion mechanisms. By contrast, in another recent paper—, 
simulations of light transmission through an array of sub- 
wavelength cylindrical holes in a tungsten layer deposited 
on glass substrate was given an interpretation based on 
dielectric guided modes. Indeed, in the wavelength do- 
main under examination, tungsten exhibits a positive 
permittivity real part**. Though SPs cannot exist, the 
transmission pattern is found very similar to that ob- 
tained with a metallic film. These theoretical results have 
recently been confirmed experimentally*'^. Surprisingly, 
the very fact that the typical transmission pattern can 
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FIG. 2: Real (ej^) and imaginary (e^) part of the chromium 
permittivity as a function of wavelength 



be observed even in non-metallic systems convinced some 
authorsi!^ to fully reject the SPs hypothesis and rather 
consider models involving non-resonant evanescent waves 
diffraction. To our knowledge, this point of view is not 
supported by recent resultsiiiiiii and the SPs must be 
replaced by other kinds of eigenmodesiiii in the non- 
metallic cases. 

The purpose of the present paper is to demonstrate 
that the transmission profile of a chromium film, in the 
restricted wavelength domain (1112 — 1292 nm) where 
the dielectric constant is positivc^^, should involve eigen- 
modes which are not SPs nor guided modes. In this case, 
SPs can be substituted by Brewster-Zennek modes (BZ 
modes) . This confirms that the observation of the effects 
described by T.W. Ebbesen et al in systems where SPs 
do not exist do not preclude the existence of a mechanism 
involving the excitation of other cigcnmodes. 



II. BREWSTER-ZENNEK MODES 

Let us recall some properties of surface modes, includ- 
ing the concept of BZ modea—^. We consider an in- 
terface between two media, 1 and 2, described by their 
respective permittivities ei and S2- It is well know that 
for an isotropic material, only p polarized modes can oc- 
cur as a surface modes. We define k as the wave vector 
component parallel to the interface and we set ko — uj/c, 
where lu is the angular frequency of the mode. Then, 
for the surface mode, the normal component of the wave 
vector is given in each medium by 



fc2 



(1) 



The subscript i denotes one of the media 1 or 2. 

In the following we enforce a real positive value of ei, 
which means that medium 1 is a lossless dielectric ma- 
terial, while we admit the form £2 = £2 + (^2 — 0) 
for medium 2, which means that this medium is allowed 
to generate losses. Applying electromagnetic boundary 
conditions pertinent to complex-response discontinuities, 
we are lead to the following well-known expression which 
controls the appearance of surface modesSi 



ei(32 + £2/?i = 



(2) 



From this, and eq. we can develop the following ex- 
pression, out of which a dispersion relation can be drawn 



j^2 £l£2 



(3) 



However, it is important to realize that in this expression, 
the factors Pi have been squared, so that eq. © and Q 
are not equivalent. The solutions of eq. (|2Jl can be found 



among those of eq. (0), but some of the solutions of eq. 
Q are spurious and should be discarded. One way to do 
this is to take the possible values of fc^ from eq. (jSjl, inject 
them in eq. and check for the correct matching. This 
actually generates a second condition, which can be used 
together with eq. Q, to make precise the acceptable 
eigenmodes branches. This condition turns out to be 



l + z 



-1 



= 
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where, somewhat surprisingly, the only variable that 
needs to be examined is the ratio z — e-il £\- The function 
of z, at the left-hand side of this equation, when carefully 
calculated, defines a border line which splits the complex 
plane into two extended connected regions. On one side 
of this border, the left-hand side of eq. (0J cannot vanish, 
while on the other side, it will vanish everywhere. The 
surface states can only show up when the ratio z lies in 
this latter region. 

The border separating the two regions can be obtained 
as follows. Both terms added in eq. ^are complex num- 
bers with the same modulus. In one of the regions men- 
tioned above, they are opposite complex numbers and 
they cancel each other; in the rest of the complex plane, 
they turn out to be identical and they cannot add to zero. 
The borderline will then be found at the complex points 
z where the sign determination of one of the square roots 
in eq. ^changes, and this requires to look for the cut of 
the complex square-root function \fZ. It is convenient to 
take this cut as the straight segment Z = —i'j (7 real and 
positive) and this choice lead us to identify two distinct 
implicit paths, 
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and 
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which form the border. This choice of the cut is dic- 
tated in part by the choice of a time-dependent oscillation 
written as exp(— iwt) and the requirement that vT = 1. 
Once the borderline is known, a quick check at the case 
z = \ shows that no cancellation can occur there, so that 
the active region where eigenmodes appear is, unambigu- 
ously, that region which does not contain z — 1. 

On the basis of this analysis, it is easily shown that, 
whatever £2 = £2 + is'i such that £2 < and £2 > 
surface eigenmodes exist. These modes are SPs and they 
will appear at metal/dielectric interfaces regardless of the 
metal component of the heterojunction. 

If £2 > 0, surface eigenmodes exist only when < £2 < 
£1 if the condition 



e'^ > £. 



I El 



£1 



(7) 



3 



FIG. 3: (Color online) Comparison between transmission and 
resonances. Solid line : a = 1000 nm, Dashed line : a — 1200 
nm. (a) : zeroth order transmission against wavelength. A 
and B are the main maxima of transmission, (b) : Resonant 
diffraction orders amplitudes related to the vectors (0, ±1) of 
the reciprocal lattice, at the substrate/metal interface (1) and 
vaccuum/metal interface (2), as a function of the wavelength 



is fulfilled. Obviously such modes are not SPs and it is 
well known that a dielectric/dielectric interface, without 
loss, cannot support a surface mode if Pi is real. Never- 
theless, surface modes can exist for a dielectric/dielectric 
interface if one of the media is a lossy material such that 
the imaginary part of its permittivity verify the above 
conditions. Such modes are known as Brewster-Zenneck 
modes. They verify Re{fc} < —y/ei and are necessarily 
related, by definition, to radiative surface modes (con- 
trasting SPs in metals such that iron, gold or copper). 
Fundamentally it is well known that SPs consist in a 
collective motion of electrons^i. In the context of an 
electromagnetic model, the electronic character of SPs 
is then contained in the permittivity £2. Indeed, the 
permittivity properties take their origin in the electronic 
characteristics of the mediumSi. This is true whatever 
the value of £3; positive or negative. So, despite the 
above-mentioned differences between SPs and BZ modes, 
both are related intrinsically to a collective electronic 
phenomenon. In this way, the transition from the SPs 
towards the BZ modes and vice versa occurs without 
discontinuities. Note that the present paper deals with 
metals only. In the case of ionic crystals in the far in- 
frared, it would be natural to expect for the BZ modes 
counterpart associated with phonon-polaritons instead of 
SPs. 



III. BREWSTER-ZENNEK-BASED 
TRANSMISSION 

Now we study the properties of the surface modes in 
the case of a chromium grating similar to those submit- 
ted to experiments by Thio et a& (see fig. In a large 
part of the spectrum of interest here, the real part of the 
chromium permittivity is negative (fig. as for simple 
metals below the plasma frequency. However, for wave- 
lengths in the interval from 1112 to 1292 nm, it becomes 
positive, as, specifically, we observe < £2 < Oj 4 and 
25, 54 < £2 < 27, 3. With vacuum and substrate (glass) 
permittivities equal to 1 and 2, 24 respectively, the condi- 
tions for the formation of BZ modes are verified at both 
interfaces of the chromium film. This is a peculiar situ- 
ation : with tungsten films in the dielectric domaini2ii&, 
for instance, the same conditions are not fulfilled and the 
BZ modes do not appear. 

The transmission spectrum can easily be understood 
if we first localize the BZ modes in the spectrum of a 
flat, homogeneous film. For normal incident light falling 



on a such a homogeneous film considered periodic on a 
square lattice of parameter a (the "empty" lattice case), 
surface modes resonance will appear at the wavelength 
Xij provided by the following equation^i 

where («,j) denotes the related vector g of the reciprocal 
lattice, such that g — (27r/a) (?, j). These wavelengths 
are indicative only, as they are based on a model where 
the cylindrical holes have not been accounted for. 

More detailed simulations which do not contain this 
simplification have been carried out. The calculations 
which will now be described are based on a coupled- 
modes method which combines a scattering matrix for- 
malism with a plane wave representation of the fields. 
This technique provides a computation scheme for the 
amplitude and polarization (s or p) of reflected and trans- 
mitted fields in any diffracted order. A brief account of 
the method has been presented elsewhere--'-^'-^^ so that 
we will not recall any technicality about these computa- 
tions. 

Following the experiment^, fig. shows the calcu- 
lated zeroth order transmission of a grating prepared by 
perforating a 100 nm thick chromium film deposited on 
glass with circular holes of 250 nm radius (fig. QJ. The 
two-dimensional square lattice parameter is a = 1000 nm 
(solid line) or a = 1200 nm (dashed line). The incident 
beam is normal to the grating surface and polarized so 
that the electromagnetic wave electric field is oriented 
along one of the holes nearest-neighbor directions {Oy 
axis). A and B on Fig. indicate the main transmis- 
sion maxima. Fig. represents the amplitude of the 
calculated resonant p-polarized diffraction orders (0,±1) 
for the substrate/metal interface (peak 1) and for the 
vacuum/metal interface (peak 2). Again, in Fig. Ob, the 
dashed line refers to a lattice parameter of a = 1200 nm. 
It can be noted that, as justified by the above discus- 
sion, these peaks do not correspond to the transmission 
maxima. We have shown that the spectral lineshape of 
the resonance can be interpreted as a Fano's profilei^ii^. 
What is remarkable here is that, for a lattice parameter 
of a = 1200 nm, the resonance has moved to a spectral 
region (1200 nm) where the real part of the chromium 
permittivity is positive, and that clearly SPs do not exist 
there. In this situation the resonance is related to the ap- 
pearance of a BZ mode. Fig 0] shows the location of both 
peaks (1) and (2) as a function of the grating parameter a 
(solid lines (1) and (2), respectively). We could compare 
these values with the wavelengths obtained from eq. |S| 
(dashed lines (1) and (2) respectively) for the appropriate 
diffraction orders. A good agreement is found between 
these calculated quantities. We show that as the lattice 
parameter increases, the resonance wavelength increases 
linearly as predicted by eq. |S| The hatched domain in 
fig. 01 is characterized by the presence of BZ modes. The 
SP domain lies everywhere else. One interesting observa- 
tion is that both kinds of surface modes (BZ modes and 
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FIG. 4: Position of both peaks (1) and (2) against grating 
parameters a. Solid lines (1) and (2) respectively are obtained 
from numerical computations. Dashed lines (1) and (2) are 
obtained from eq.(8). A (circle dots) and B (square dots) 
: location of both main maxima of transmission against a. 
Hatched domain : BZ modes, elsewhere : SPs modes. 

SPs modes) substitute each others in a continuous way 
in accordance with our previous explanation in section 
IL The location of transmission maxima A (circle dots) 
and B (square dots) are also shown on fig. 0] These, 
as already underlined, do not indicate directly the eigen- 
modes resonance wavelengths though these transmission 
maxima are known to be related to the surface modes via 
a Fano profile. 

Note that, in a previous worl^i, it was inferred that the 
SPs wavelength were red-shifted, compared to those cal- 
culated for an empty lattice. It was assumed that such a 
shift was needed to make the transmission maxima coin- 
cide with SPs wavelengths and the justification suggested 
for this shift was based on the presence of holes. This is 
not observed in actual detailed calculations. The present 
results show that the exact resonant wavelengths of the 
surface modes, as computed numerically, are not shifted 
enough to justify their location on transmission maxima. 
Nevertheless, as mentioned before, an eigenmode reso- 
nance does not necessarily coincide with a transmission 
maximum, but a maximum-minimum pair indicates the 
presence of a nearby eigenmode resonance^"^. In fact, in 
this case, the surface modes resonances are incidentally 
located closer to the minima of the transmission. It is 
interesting to note that Thio et al have compared exper- 
imentally the transmission from gratings similar to ours, 
with lattice parameters a = 1000 nm and a = 1200 nmS. 
The present computations of the transmission matched 



very well with these datck^. We take this as a strong evi- 
dence of the relevance of the Brewster- Zennek modes in 
the description of these experiments. 



IV. CONCLUSION 

In the context of studies of the optical transmission 
properties of subwavelength hole arrays, we have shown 
that Brewster-Zennek modes can substitude surface plas- 
mons when considering a layer made from a lossy dielec- 
tric medium, for specific lattice parameters. When mod- 
elling a grating etched in a chromium film we can use the 
concept of Brewster-Zennek modes to explain the spec- 
tra found in previous experiments. This study recalls 
the existence requirement of eigenmodes for observing 
the typical transmission lineshapes and we have shown 
that, contrasting a frequent believe, surface modes wave- 
lengths based on an empty lattice model are not significa- 
tively red-shifted to justify their location on transmission 
maxima. 
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